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Pion and photon production cross sections are analyzed in proton-proton and proton-nucleus colli- 
sions at energies 20 GeV < y/s < 60 GeV. We separate the proton-proton and nuclear contributions 
to transverse-momentum broadening and suggest a new mechanism for the nuclear enhancement in 
the high transverse-momentum region. 

PACS number(s): 12.38.Bx, 13.85.Ni, 13.85.Qk, 24.85. +p, 25.75.-q 



(N 
> 

o 

m 
O 
as 

OS 



o 

=5 



X 



Experimental results on pion and photon production 
in high-energy hadron-nucleus collisions show an extra 
increase at high transverse momentum (j>t) over what 
would be expected based on a simple scaling of the ap- 
propriate proton-proton (pp) cross sections. The nuclear 
enhancement is referred to as the Cronin effect , and is 
most relevant at moderate transverse momenta (3 GeV 
^ Pt ^ 6 GeV) 0] . In relativistic nuclear collisions this 
momentum region is at the upper edge of the px win- 
dow in Super Proton Synchrotron (SPS) experiments at 
y/s = 20 AGeV, and is measurable at the Relativistic 
Heavy Ion Collider (RHIC) at y/s = 200 AGeV. The 
importance of a better theoretical understanding of the 
Cronin effect continues to increase as new data ap- 
pear from the SPS heavy-ion program and as the commis- 
sioning of RHIC approaches. This calls for a systematic 
study of particle production moving from pp to proton- 
nucleus (pA) collisions. In the present work we analyze 
the Cronin effect in inclusive 7r° and 7 production. 

In the past two decades the perturbative QCD (pQCD) 
improved parton model has become the description of 
choice for hadronic collisions at large px- The pQCD 
treatment of hadronic collisions is based on the assump- 
tion that the composite structure of hadrons is revealed 
at high energies, and the parton constituents become the 
appropriate degrees of freedom for the description of the 
interaction at these energies. The partonic cross sections 
are calculable in pQCD at high energy to leading order 
(LO) or next-to-leading-order (NLO) §-|). The parton 
distribution function (PDF) and the parton fragmenta- 
tion function (FF), however, require the knowledge of 
non-perturbative QCD and are not calculable directly 
by present techniques. The PDFs and FFs, which are 
believed to be universal, are fitted to reproduce the data 
obtained in different reactions. In recent NLO calcula- 
tions the various scales (Q,Aqcd, etc.) are optimized to 
improve the agreement between data and theory PJlC|]. 

In theoretical investigations of 7r° and 7 production in 
pA collisions another method appeared and became pop- 



ular |1 1|-|13[ : the different scales of the pQCD calculations 
are fixed and the NLO pQCD theory is supplemented by 
an additional non-perturbative parameter, the intrinsic 
transverse momentum (kx) of the partons. The presence 
of an intrinsic transverse momentum, as a Gaussian type 
broadening of the transverse momentum distribution of 
the initial state partons in colliding hadrons was inves- 
tigated as soon as pQCD calculations were applied to 
reproduce large-pr hadron production [jl4 15 1. The av- 
erage intrinsic transverse momentum needed was small, 
(kx) ~ 0.3 — 0.4 GeV, and could be easily understood in 
terms of the Heisenberg uncertainty relation for partons 
inside the proton. This simple physical interpretation 
was ruled out as the only source of intrinsic kx by the 
analysis of new experiments on direct photon production, 
where (kx) ~ 1 GeV was obtained in the fix target Teva- 
tron experiments p] p"3| and (kx) ~ 4 GeV was found at 
the Tevatron collider for muon, photon and jet produc- 
tion [|l3|. New theoretical efforts were ignited to under- 
stand the physical origin of (kx) Parallel to these 
developments, kx smearing was applied successfully to 
describe ultrarelativistic nucleus-nucleus collisions |}|4|] 
and J/ip production at Tevatron and HERA One 
possible explanation of the enhanced &T-broadening is 
in terms of multiple gluon radiation [jl7| , which makes 
(kx) reaction and energy dependent. In the absence of 
a full theoretical description, intrinsic kx can be used 
phenomenologically in pp collisions. A reasonable repro- 
duction of the pp data is a prerequisite for the isolation 
of the nuclear enhancement we intend to focus on. 

In the lowest-order pQCD parton model, direct pion 
production can be described in pp collisions by 

E -K~jt~ = ^2 J dx ~L dx 2 fa/p(xi,Q 2 ) fb/p(x2,Q 2 ) 
abed 

K^(ab^cd) D */°M 2 ) , (l) 

dt 7TZ C 

where f a /p(x,Q 2 ) and fb/p(x,Q 2 ) are the PDFs for the 
colliding partons a and b in the interacting protons as 



1 



functions of momentum fraction x and momentum trans- 
fer Q, and a is the LO hard scattering cross section 
of the appropriate partonic subprocess. The K-factor 
accounts for higher order corrections |fl5|] . Comparing 
LO and NLO calculations one can obtain a constant 
value, K w 2, as a good approximation of the higher 
order contributions in the px region of interest |Q. In 
eq.(|l]) D 7r / C (z c , Q 2 ) is the FF for the pion, with the scale 

Q = Pt I ' z Cl where z c indicates the momentum fraction of 
the final hadron. We use a NLO parameterization of the 
FFs |2(|. Direct 7 production is described similarly 0. 

The generalization to incorporate the kx degree of free- 
dom is straightforward HQ. Each integral over the par- 
ton distribution functions is extended to /c^-space, 

dx f a /p(x, Q 2 ) -> dx d 2 k T g(k T ) f a / p (x, Q 2 ) , (2) 

and, as an approximation, g(kr) is taken to be a Gaus- 
sian: g{kr) — exp(— k T /(k T ))/n(k T ). Here (k T ) is the 
2-dimensional width of the kx distribution and it is re- 
lated to the average transverse momentum of one parton 
as (k 2 ) = 4(fc T ) 2 /7r. 



(k T ). The error bars display a A 2 = A 2 ^ ± 0.1 un- 
certainty in the fit procedure. The uncertainty is small 
at y/s = 20 — 30 GeV and relatively large at y/s « 60 
GeV, indicating that (k T ) is more sharply determined at 
lower energies. We use a value of (k T ) = 3 GeV 2 for tt° 
at y/s = 27 A GeV. The value of (k T ) appears to increase 
with energy. The dashed line serves to guide the eye and 
indicates a linear increase to a value of (kr) = 3.5 GeV 
at y/s = 1800 GeV @. 

Furthermore, we analyzed the data from pp — > jX 
reactions [ p5| , p6| , p9| ^33]| . These results are also shown 
in Fig.l. In this case much lower values are obtained 
for (k T ) with large uncertainty. The dotted line repre- 
sents the results obtained in Ref. Jll| from diphoton and 
dimuon data. In the following we use (k T ) =1.2 and 1.5 
GeV 2 at y/s = 31.6 and 38.8 GeV, respectively. 

Fig. 2. displays the ratio of data to our calculations as 
a function of xt — 2px/-\/s for pp — > tt°X, applying the 
obtained best fit values for (k T ) in the different experi- 
ments. For pp — ► 7V data we achieve a similarly good 
agreement. Further details of the comparisons with data 
and with other calculations will be published elsewhere 
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FIG. 1. The best fit values oUtft) in pp -> n°X ^^-|| 
and pp ■yX reactions [ p5|^ , pi| -p^] (upper panel), and a 
calculated (&t) (lower panel). See text for the lines. 
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FIG. 2. The ratio of data to theory for pp — > ty°X reac- 
tions Jl||2^-||l applying the best fit for (fe|), x T = 2p T /y/s. 
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We applied this model to describe the measured data 
tt°X reactions |l21-2(|. If 7r + and n~ pro- 
duction was measured, we constructed the combination 
7T° = (tt+ + 7T~)/2, given by the FFs |o| we use. The 
calculations were corrected for the finite rapidity win- 
dows of the data. The Monte-Carlo integrals were carried 
out by the standard VEGAS-routine |27|. For the PDFs 



we used the MRST98 set 28 , which incorporates an in- 
trinsic kr- The scales are fixed, Kjjglrif = 4) = 300 
MeV and Q — pt/2. We fitted the data minimizing 
A 2 = y^/Daba — Theory) 2 /Theory 2 in the midpoints 
of the data. Fig.l. shows the obtained fit values for 



With 7T° and 7 production in pp collision reasonably 
under control, we turn to the nuclear enhancement in pA 
collisions. In minimum-biased pA collisions the pQCD 
description of the inclusive pion cross section is based on 

E„ ^ =2J Jd 2 bt A (b) J dx 1 dx 2 d 2 k Tta d 2 kT,b g(kr,a) 



'ibed 



g(k T ,b)fa/p(x 1 ,Q 2 )f b / p (x 2l Q 2 ) K 



da D v/c (z c ,Q 2 ) 



dt. 



(3) 



Here b is the impact parameter and ^(6) is the nuclear 
thickness function normalized as J d 2 b tAib) = A. For 
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simplicity, we use a sharp sphere nucleus with t A (b) = 
2p ^n{- b 2 , where R A = 1.14A x /3 and p = 0.16 
fm~ 3 . Next we discuss the nuclear enhancement of (k 2 -). 

The standard physical explanation of the Cronin ef- 
fect H is that the proton traveling through the nucleus 
gains extra transverse momentum due to random soft 
collisions and the partons enter the final hard process 
with this extra kx . In our approximation initial soft pro- 
cesses increase the value of but this effect does not 
depend on the scale, Q 2 , of the hard process to occur 
later. However, it may depend on the initial center-of- 
mass energy, y/s. Furthermore, it is important to note, 
that in our description not all participant protons are au- 
tomatically endowed with the extra enhancement, 
only the parton distribution of the colliding protons is 
affected according to the number of soft collisions suf- 
fered. To characterize the (k 2 -) enhancement, we write 
the width of the transverse momentum distribution of 
the partons in the incoming proton as 

(J4) P A = {kr)p P + C ■ h pA (b) . (4) 

Here h p A{b) describes the number of effective nucleon- 
nucleon (NN) collisions at impact parameter b which im- 
part an average transverse momentum squared C . 

Naively all possible soft interactions are included, 
hp l X(b) = VA{b) - 1, where v A (b) = (Jnn t A (b) is the 
collision number at impact parameter b with <jnn be- 
ing the total inelastic cross section. Applying this model 
to 7T° production in pA (A = Be,Ti,W) collisions at 
y/s = 27 A GeV 0, we extract C$ e = 0.8 ± 0.2 GeV 2 , 
i = 0.4 ± 0.2 GeV 2 , and = 0.3 ± 0.2 GeV 2 . The 
target dependence of the extra enhancement in (k 2 -,) is in- 
consistent with the assumption of a target-independent 
average transverse momentum transfer per NN collision 
and/or with the form of h^(b). 

Inspired by the A-dependence, we propose another 
physical picture of the nuclear enhancement effect. Ac- 
cording to this mechanism, the incoming nucleon first 
participates in a semi-hard (Q 2 ~ 1 GeV 2 ) collision re- 
sulting in an increase of the width of its kr distribu- 
tion. There is at most one collision able to impart this 
increased value of (fc 2 ,), characteristic of the no-longer 
coherent nucleon. Further soft or semi-hard collisions do 
not affect the kx distribution of the partons in the incom- 
ing nucleon. This saturation effect can be approximated 
well by a smoothed step function (b) defined as 

f if v A {b) < 1 

h$(P) = { - 1 if 1 < u A (b) < 2 . 

[ 1 if 2 < v A (b) 

The saturated Cronin factor is denoted by C sat . 

Using the same pA data as previously Q, C sat = 1.2 
GeV 2 gives a good fit for all three targets, Be, Ti and 
W. Fig. 3. displays our results with (full line) and with- 
out (dashed line) the Cronin enhancement. The lower 



panel shows the data/theory ratio on a linear scale for 
the pW case. It is interesting to note that with the sat- 
urated Cronin effect the remaining deviation from one 
in the data/theory ratio is similar in shape to the pr- 
dependent K factor obtained in Ref. [0]. The pBe — ► n°X 
data at energies y/s = 31.6 GeV and 38.8 GeV Q are 
also described quite well with the same saturating Cronin 
effect §§. 
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FIG. 3. Cross section per nucleon in pA — > 7r X reactions 
(A = Be,Ti,W), data from §. We show C sat = 1.2 GeV 2 
(full lines) and C sat = (dashed lines). The lower panel 
displays the data to theory ratio on a linear scale for the pW 
collision. 

Let us now discuss photon production in the same en- 
ergy region. We speculate that the nuclear enhancement 
does not depend on the outgoing particle, and thus use 
the same C sat for 7 as for ir° production. Fig. 4. shows 
the pBe — > 7V reaction at energies y/s = 31.6 GeV and 
V* = 38.8 GeV. The data are from Ref. @. The calcu- 
lations are carried out with C sat = 1.2 GeV 2 (full lines). 
For comparison we show the results without nuclear en- 
hancement, C sat =0 (dashed lines). In the lower panel 
the data/theory ratio is displayed for y/s = 31.6 GeV. 

The common value of C sat in the studied energy range 
indicates that the extra enhancement in pA collision is 
independent of the produced final state particle. We 
interpret this C sat as the square of the characteristic 
transverse momentum imparted in one semi-hard colli- 
sion prior to the hard scattering. Independent of the 
details of the mechanism, the extra nuclear enhancement 
in eq. (4) appears to have the same total value, which is 
on the scale of the intuitive dividing line between hard 
and soft physics, C ■ h pA (Q) w 1 GeV 2 . For RHIC predic- 
tions it would be necessary to see the energy dependence 
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of C sat , and whether the same nuclear enhancement is 
obtained at much higher energies. 




Px (GeV) 

FIG. 4. Cross section per nucleon in the pBe —* 7A 
reaction at y/s = 38.8 GeV (dots) and at y/s — 31.6 GeV 
(full triangles) [jfj. Solid lines indicate C 3at = 1.2 GeV 2 , 
dashed lines mean C sat = 0. Lower panel shows data/theory 
for yfl = 31.6 GeV. 

In the present letter we separated the pp and nuclear 
contributions to the width of the parton transverse mo- 
mentum distribution, and proposed a saturating model 
of the Cronin effect in proton- nucleus collisions. Accord- 
ing to our picture, the incoming proton suffers at most 
one semi-hard scattering prior to the hard parton scatter- 
ing. In the semi-hard collision the width of the transverse 
momentum distribution of the partons inside the incom- 
ing proton increases. This prescription describes the px 
dependence of the nuclear enhancement in tt° and 7 pro- 
duction remarkably well. We are confident that complete 
NLO calculations, which are left for future work, will 
improve the agreement further. Systematic pA experi- 
ments are needed to determine the energy dependence of 
C sat (s) and to extrapolate to AA collisions at RHIC. 
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